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Finite-size scaling analysis turns out to be a powerful tool to calculate the phase diagram as well 
as the critical properties of two dimensional classical statistical mechanics models and quantum 
Hamiltonians in one dimension. The most used method to locate quantum critical points is the so 
called crossing method, where the estimates are obtained by comparing the mass gaps of two distinct 
lattice sizes. The success of this method is due to its simplicity and the ability to provide accurate 
results even considering relatively small lattice sizes. In this paper, we introduce an estimator 
that locates quantum critical points by exploring the known distinct behavior of the entanglement 
entropy in critical and non critical systems. As a benchmark test, we use this new estimator to 
locate the critical point of the quantum Ising chain and the critical line of the spin-1 Blume-Capel 
quantum chain. The tricritical point of this last model is also obtained. Comparison with the 
standard crossing method is also presented. The method we propose is simple to implement in 
practice, particularly in density matrix renormalization group calculations, and provides us, like the 
crossing method, amazingly accurate results for quite small lattice sizes. Our applications show that 
the proposed method has several advantages, as compared with the standard crossing method, and 
we believe it will become popular in future numerical studies. 
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I. INTRODUCTION 

In recent years, the connection between the quantum 
correlations and the entanglement properties of quan- 
tum critical systems motivated the research on the pos- 
sible characterization of critical phenomena via quantum 
information concepts i - — In particular, several of these 
studies have been devoted to the identification of possi- 
ble measures of entanglement that give a precise local- 
ization of the quantum critical point in quantum chains. 
These studies basically considered local measures of en- 
tanglement like the concurrence and the negativity^ - — 
or other measures of quantum correlations that are not 
always simple to calculate*^— The concurrence and the 
negativity are easy to measure, however can not be con- 
sidered universal measures, since depending on the dis- 
tance of the sites they are measured, and the model, they 
may detect or not the critical points^ ' 2 ' 9 ' 17 

On the other hand, the standard theory of finite-size 
scaling (FSS) gives us a quite simple way to locate quan- 
tum phase transitions. The procedure, usually named as 
the crossing method (CM), is based on the crossing of 
the functions LG(L) = L'G(L') at the critical points^ 
where G(L) and G(L') are the mass gaps of the quantum 
chains with lattice sizes L and L' , respectively. This is 
the most used method to locate a phase transition in a 
quantum chain. The success of this method is due to the 
fact that it provides precise estimates of critical points 
even considering small system sizes. This point is very 
important because, with the exception of the exactly sol- 
uble models, we do not have access of the ground state 



properties of large systems. Normally, for a given lattice 
size, the CM method gives the most accurate localization 
of critical points. 

A natural question concerns the existence of an uni- 
versal estimator, based on the changing of entanglement 
at the quantum critical points, that gives the same order 
of precision, for a given numerical effort, as the standard 
CM. In this paper, we are going to present such estimator 
for the quantum chains. The estimator will be obtained 
by exploiting the well known violation of the area law of 
the entanglement entropy of a subsystem with the rest 
of the quantum chain that happens in a critical quantum 
point. 

The entropy itself could, in principle, be used to locate 
the quantum critical points. However, as we are going 
to show in two examples, we need to use much larger 
lattices, to achieve the same precision, as compared with 
the CM. The estimator we propose is given by the dif- 
ference of the entanglement entropy when we split the 
quantum chain in subsystems of different sizes. As a 
benchmark test, we are going to use this estimator to 
determine the localization of critical points of two well 
known models: the quantum Ising chain and the spin- 
1 Blume-Capel model (BCM). These are good examples 
since in the Ising model we have a well known critical 
point, and in the BCM we possess a line of second order 
phase transition in the Ising universality, that ends in a 
tricritical point in the universality class of the tricritical 
Ising model. Surprisingly, the estimates for the critical 
and tricritical points obtained by using the proposed es- 
timator are as precise as the standard CM, for the same 
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lattice sizes. Moreover, the new method in comparison 
to the CM has some advantages. We need only to calcu- 
late the ground state eigenf unction, instead of the mass 
gaps as in the CM. Furthermore at the same time as the 
quantum critical point is located the conformal anomaly 
of the underlying conformal field theory governing the 
critical point is also evaluated. Our applications also 
shows that in the case of first order transitions (gapped) 
the proposed estimator, differently from the CM, shows 
a distinct finite-size behavior already for quite small lat- 
tices, that in principle can be used to detect if the phase 
transition is discontinuous (first order) or not. 

The paper is organized as follows. Since the proposed 
method to locate quantum phase transitions is going to 
be compared with the standard CM, we present both 
methods in the next section. In the sections III and IV, 
we present our results for the quantum Ising chain and 
the BCM, respectively. Finally, in section V we summa- 
rize our results and present our conclusions. 



II. ESTIMATORS OF QUANTUM CRITICAL 
POINTS 



Most of the critical chains besides being scale invariant 
are also conformal invariant, its long-distance physics is 
governed by a conformal field theory (CFT) with cen- 
tral charge c. In the scaling regime 1 << £ << L the 
entanglement entropy Sl {£, A) of the ground state be- 
haves differently if the system is critical (A = A c ) or not. 
It increases logarithmically with the correlation length 
Cl(A), or equivalently, with the decrease of the mass gap 
G L (A) = 1/Cl (A) , namely^^-^ 



S L (£,X) 



' 7 ln[|sin(f)] +/3, A = A C 
- 7 ln[G L (A)]+/3', A^A C 



(2) 



where 7 is an universal constant related with the central 
charge c, namely 7 = c/3 for periodic chains and 7 = c/6 
for chains with open boundaries. The constants /3 and 
/3' are non universal and model dependent. Here, we 
consider only PBC in order to avoid boundary effects. 

The estimator we propose to locate a quantum criti- 
cal point A c is given by the coupling A^ that gives the 
maximum value of the entanglement entropy difference 
(MVEED) given by: 



Let us consider a quantum chain whose Hamiltonian 
H(X) has a quantum critical point at A = A c . The FSS 
theory tell us that at the quantum critical point, the finite 
chain with L sites will have a mass gap Gz,(A) ~ as 
L — > 00. This gives us the most popular method used 
to locate phase transitions. The finite-size estimates of 
the critical point {\l,L') are obtained from the mass gaps 
of two distinct lattice sizes. They are obtained from the 
crossing of the function Fl{X) = LGl{X) for two lattice 



sizesi 1 ^ 1. e. 



LG l (X)\ XlM =L'G v (X)\ Xi 



(1) 



For this reason, this method is also known as the CM. 
Its application is quite simple because it is necessary to 
evaluate only the two lowest eigenvalues of H(\). A diffi- 
culty of this method happens in the case where the quan- 
tum chain undergoes a first-order phase transition with 
a small (but finite) mass gap Goo. In this case, we also 
find crossing points Xl,L', since for L < I/Goo we should 
expect Gl(A) ~ l/L, as in a true critical point. 

The purpose in this paper is to introduce and test a 
alternative method to locate critical points in quantum 
chains. The method we propose exploits the changing of 
behavior of the entanglement entropy of a segment of the 
quantum chain, as we cross a quantum critical point. 

In order to present our method, let us first define the 
well known entanglement entropy. Consider a quantum 
chain with L sites, described by a pure state whose den- 
sity operator is p. Let us consider that the system is 
composed by the subsystems A with £ sites (£ — 1, . . . , L) 
and B with L — I sites. The entanglement entropy is de- 
fined as the von Neumann entropy Sl (£) = — Trp.4 In p^, 
associated to the reduced density matrix pj^ = Trgp. 



AS L (X) = S L (£,X)-S L (£',X), 



(3) 



with £ = L/2 and £' = L/A. According to Eq. as 
L — > 00, we should have 



A c 

A/ A c 



(4) 



This means that as L increases the estimator goes to 
zero, except at the critical point, where its value is pro- 
portional to the central charge c. Although in the defini- 
tion it is assumed that the lattice sizes are multiples 
of 4, for other lattice sizes a simple variation of Eq. ([3]) 
can also be defined. In principle, we could choose two ar- 
bitrary subsystem sizes £ and £' . However, the signature 
of the violation of the entropic area law is observed more 
clear for £,£' >> 1. Then, we have to choose £,£' ~ L. 
On the other hand, we also want to distinguish ASl (£, £') 
from zero. This means that we have to choose small val- 
ues of £' , when compared with I ~ L, and satisfying 
£' >> 1. So, it seems to be natural to choose £ = L/2 
and £' = L/4. 

In the following sections, we are going to use the 
MVEED estimator to locate the critical points of two 
different models, and compare them with those obtained 
through the standard CM. Note that in the same line, we 
could, in principle, also use as estimators the difference 
of the a-Renyi entropies. However, for a > I unusual 
corrections appears^ and the finite-size effects became 
more relevant (see also Refs. [24ll25| ). 

The key ingredient in the MVEED method is the vio- 
lation of the entropic area law at the critical point. This 
violation have been observed in several systems such as 
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the conformal invariant critical chains. We could naively 
expect that the entanglement entropy of a c£-dimensional 
system would behave as S ~ L d , since the entropy is 
usually an extensive quantity. However this is not cor- 
rect, the entanglement entropy of the subsystems A and 
B, with arbitrary volumes, are identical, i. e.: Sa — 
-TrpA^ Pa = Sb = — Trpglnpg. Since both subsystems 
share the same area, we should expect the area-law de- 
pendence of the entropy, i. e., Sa — Sb ~ £,d-i 2^-28 
This means that the information is shared only among 
the degrees of freedom localized around the surface sep- 
arating both systems. However, at the critical point the 
correlation length diverges and both subsystems share 
much more information among themselves. We thus ex- 
pect, in this case, a violation of the entropic area law. As 
we already mention, the one dimension conformal invari- 
ant critical systems [see Eq. (2)] are celebrated exam- 
ples where this law is violated. There are also examples, 
in dimension d >1, where the entropic area law is vi- 
olated, such as some gapless fermionic systems with a 
finite Fermi surfaced - — It is interesting also to point 
out that in some disordered one-dimensional systems the 
entropic law can be violated in a different way than the 
one presented in Eq. ([2])<2& 

We should mention that the difference of the entangle- 
ment entropy of chains with distinct sizes were used to 
extract the central charge c in one dimension quantum 
systems with good accuracy^ - — Lauchli and Kollath in 
Ref. Hi] also used a similar procedure to detect the contin- 
uum line of critical points of the bosonic Hubbard chain. 
In their procedure they considered the following estima- 
tor: AS LK (L, V = L/2, A) = S L (L/2, A) - S L/2 (L/4, A). 

For the above model, and also for the XX Z quantum 
chain, where a continuum line of critical points appears 
by tuning a single coupling, the difference of entropies 
tends toward a constant, appearing a plateau along the 
critical line. In the non critical region they go to zero, as 
the lattice size increases. In this case, we should not use 
the maximum of ASl(A) [or the maximum of ASlk (A)] 
but its over all size dependence to locate the endpoint 
of the critical line. We also have calculated ASl(A) for 
the XX Z chain, for L < 100, and we confirmed that the 
end point of the critical line is better estimated by the 
overall sizes dependence of ASl(A), as compared with 
the MVEED. 

For the sake of comparison, we also test the estima- 
tor proposed by Lauchli and Kollath to detect the criti- 
cal point of the quantum Ising chain. As our estimator, 
ASlk(L, V , A) also has a maximum close to the criti- 
cal point [see next section and Fig. 1(e)]. So, it is also 
possible to estimate the critical point by considering the 
coupling that gives the maximum values of ASlk(L, L'). 
However, note that to obtain one estimate of the critical 
point, with this estimator, we need to consider two dis- 
tinct lattice sizes, one being twice the size of the other. 
This demand larger numerical effort. Besides that we 
verify, for both models with L = 16, that the errors of 
our estimates are approximately two orders of magnitude 



A 75 
5 



0.5 



(a) 






f* 2 ^— X=1.0 




— X=1.01 








L=64 







0.8 



10 



20 



30 



(b) 


X=1.0 




■ - x=l . 1 








L=12 




0.4 
0.7 



> 3 


4 e 5 6 






*** Ising Model 


i — ^> 

- . 




— L=8 




■ — L=16 


lV''i„ liii 


■ -L=12 




«L=24 





0.9 



1.1 



0.06 




Figure 1: (Color online), (a) The entanglement entropy 
Sl{£, A) vs £ for the Ising model for a system with size L = 64 
and some values of A (see legend), (b) Same as (a) but for 
L — 12. (c) Sl(L/2, A) vs A for several values of L. (d) 
dS L (L/2,\)/d\ vs A. (e) ASl(A) vs A. 



smaller than the ones obtained by considering the esti- 
mator AS LK (L, L' = L/2, \)M 



III. QUANTUM ISING CHAIN - CRITICAL 
POINT DETERMINATION 



The Ising quantum chain describes the dynamics of 
spin- 1/2 localized spins whose Hamiltonian is given by 



H 



Ising 



(5) 



where a x , a z are spin-1/2 Pauli matrices. It depends 
on the parameter A, and for simplicity, hereafter we are 
going to consider only periodic chains. This model has 
a critical point A c = 1 that can be obtained from its 
exact solution, or even more simply from its self dual 
property^ 

Let us show first the difficulty of extracting the lo- 
calization of the critical point by directly using the Eq. 
((2|). In Fig. 1(a), we present for the lattice size L = 64 
the entanglement entropy Sl(£, A) of the Ising model for 
A = 1.0, 1.01 and A = 1.1. These results can be obtained 
from free fermion technique i 44 i 45 We also use these results 
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to test the precision of our density matrix renormaliza- 
tion group (DMRG) calculation. For the critical coupling 
A = A c = 1 the entanglement entropy has a quite good 
fit with the expected behavior given in Eq. ([2]) with the 
values c = 0.5003 and (3 — 0.4781, which are close to the 
exact ones4& On the other hand, for a point clearly away 
from the critical point, like A = 1.1, the entropy tends 
towards a constant as I increases, as expected from Eq. 
([2]). Although we see a clear distinct behavior of the 
entropy for such anisotropies, for anisotropies closer to 
the critical point it is quite difficult to distinguish such 
behavior even for relatively large lattices. For example, 
as we can see in Fig. 1(a), it is difficult to discern if we 
are at the critical point or not at A = 1.01 for a chain 
with lattice size L = 64. In fact, for this coupling, we 
also obtain a nice fit with the predicted critical behavior 
given in Eq. ([2]). The fit to our numerical data gives us 
c = 0.489 and (3 = 0.472, which are close to the values 
at the critical point. The reason, for not seeing a sat- 
uration, as we increase £, for A = 1.01, is that at this 
coupling and for such lattice size the correlation length 
Ci(A) = 1/Gl(A) > L. Actually, it is possible to dis- 
tinguish the critical and the non critical behavior only 
when one consider lattice sizes L >> C,. This means that 
for relative small lattices (L ~ 20), differently from the 
CM, it is not possible to obtain reasonable estimates for 
the critical coupling constants only from the behavior of 
Sl(£, A). Indeed, as shown in Fig. 1(b), for the lattice 
size L = 12, it is even difficult to distinguish the behav- 
ior at the anisotropies A = 1 and A = 1.1. Certainly, the 
distinct behavior of Sl{£, A) as a function of £, should be 
more evident at I = L/2. So, let us check if the signa- 
ture of the critical point is more evident for I — L/2. In 
Fig. 1(c), we present Sl(L/2, A) as a function of A for lat- 
tice sizes L = 8, 12, 16, 24 and 36. In this figure only for 
L > 24 we start to see a clear change of behavior around 
the expected critical point. We observed that the maxi- 
mum occurs at A ~ 0.97. This means that the simple use 
of Sl{L/2,X) will give an acceptable estimate only for 
quite larger lattices, as compared with those necessary 
to obtain reasonable estimates through the use of the 
CM. The fact that we need to consider large systems to 
observe in Sl(L/2,A) the signature of the critical point 
is not particular of the Ising model. This behavior have 
been observed also in other model a 39 ' 47 i 4s [see also Fig. 
2(a)]. We should also mention that we did not observe 
a clear signature of the quantum critical point, for small 
lattice sizes, in the derivative of the entanglement entropy 

dS(L/2,X) r -p. 

dx ■ [see Fig. 1(d)]. 

Now, we are going to show that ASl (A) , defined in Eq. 
already presents the signature of the phase transition 
for small system sizes. This fact is highly desirable, be- 
cause it will not need a huge numerical effort to detect 
the localization of the phase transition. In Fig. 1(e), we 
depict the difference ASl(A) as a function of A, for the 
same lattice sizes presented in Fig. 1 (c) . We notice clearly 
two interesting features: (i) the maximum of ASl (A) ap- 
pears close to the true critical point A c = 1, and more 



L c L 
8 0.99307 0.50374 
(1.00197) 

12 0.99813 0.50137 
(1.00049) 

16 0.99939 0.50073 
(1.00014) 

24 0.99984 0.50031 
(1.00005) 

36 0.99996 0.50014 
(1.00001) 

48 0.99997 0.50008 
(1.00001) 

exact 1.0 0.5 

Table I: The finite-size estimates of and c L for the quantum 
Ising model acquired from the MVEED method and CM. The 
results in parentheses are from the CM. The exact values of 
these quantities are also shown. 

important (ii) the signature of the phase transition is ob- 
served already for small lattices, as happens in the CM 
(see Table I). 

In Table I, we present the finite-size estimates of A^ ob- 
tained from the maximum value of ASl (A) . We also show 
in this table the finite-size estimate of the central charge 
c L = 6ASl(A^ )/ln(2). As we can note, it is grateful to 
see that accurate estimates, which are around 0.01% of 
the exact value (A c = 1, c = 1/2), are already obtained 
for the small lattice size L = 8. We also show in Table I, 
for comparison, the estimated values acquired from the 
standard CM [see Eq. (JTJ]. We learn from this applica- 
tion that the results obtained with the MVEED method 
are as precise as the standard CM. Differently from the 
CM, our method has the advantage to give us, as a bonus, 
good estimates of the conformal anomaly c at the same 
time we got the estimate \\. In order to see that the 
success of this method in not a particular case of the ex- 
actly integrable quantum Ising chain, in the next section 
we present results for a non integrable quantum chain. 



IV. BLUME-CAPEL QUANTUM CHAIN - 
CRITICAL POINT DETERMINATION 

The Blume-Capel quantum chain is obtained by the 
time-continuum limit of the Blume-Capel model in two 
dimensions. It describes the dynamics of spin-1 localized 
particles, with Hamiltonian given by 

HBC = -J2(s*s* +1 -6(s*f- 7 s*), (6) 

3 
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where s x and s z are the spin-l operators. This Hamilto- 
nian, differently from the quantum Ising chain, has two 
coupling constants (6 and 7) and is not exactly inte- 
grable. The phase diagram in the plane S — 7 is known 
from earlier numerical studies based in the CM (see Fig. 
1 of Rcf. 49). For values 7 > 7t r the Hamiltonian 
has a quantum critical line <5 C (7) governed by a CFT in 
the same universality class of the quantum Ising chain 
(c = 1/2). At 7t r the model has a quantum tricritical 
point at S tr in the universality class of the tricritical Ising 
model, having central charge c = 7/10. For 7 < 7 tr there 
is a line ^(7) of first-order phase transitions (gapped). 
The numerical estimate of the tricritical point (7t r ,^tr) 
was acquired by a generalization of the CM. It was ob- 
tained either by the crossing of the two lowest mass gaps 
of two lattice sizes, or by the crossing of the first mass 
gap of three lattice sizes. These previous results were 
obtained for small lattice sizes (L < 9) and the mass 
gaps were calculated by using the Lanczos method. 50 In 
this work, we also used Lanczos method together with the 
CM to extend the previous finite-size estimates of (j c , S c ) 
up to lattice size L = 15. Some of these estimates are 
presented in Table II. 

Now, let us calculate the critical line of the model by 
using the MVEED [Eq. (gj)]. Although it is also possible 
to calculate the entanglement entropy with the Lanczos 
method, here, we used the DMRG technique 5 -^ to do this 
task. In the DMRG the entanglement entropy is easily 
obtained since the reduced density matrix is naturally 
calculated in all sweeps of the DMRG. In our DMRG, we 
kept up to m = 600 states per block in the final sweep. 
We have done ~ 4 — 6 sweeps, and the discarded weight 
was typically 10~ 12 — 10~ 15 at that final sweep. 

Like in the case of the quantum Ising chain, let us con- 
sider first the behavior of the entanglement entropy for 
the subsystem with size 1 = L/2, i. e., Sl(L/2, 7, 5). We 
consider initially a value of 7 where a second-order phase 
transition in the Ising universality is expected. This is 
the case for 7 = 1.1, where a critical point happens for 
S c ~ 0.3135 (see Table II ). 

In Fig. 2(a), we show the entropy Sl(L/2, 7 = 1.1, S) 
as a function of 5. Similarly as in the application for 
the quantum Ising chain, this entropy does not present 
any trace of a quantum phase transition for small system 
sizes. We start to observe the signature of the phase tran- 
sition only for relative large lattice sizes (L > 24) . Similar 

results is also seen in the derivative dSL ^] ,S ^ , as depicted 
in Fig. 2(b) for 7 = 1.1. On the other hand, if we now 
use our estimator AS/, (7, S), as shown in Fig. 2(c), we al- 
ready see a clear appearance of a peak for the small lattice 
size L = 8. The finite-size estimates <5£, for three values 
of 7, obtained from the MVEED method are presented 
in Table II. It is also given in this table, for comparison, 
the values of S' 



L-2,L 



obtained through the standard CM 
[Eq. (p])]. We clearly see in the table that the results 
acquired from the MVEED method converged already to 
almost five digits for small lattice sizes like L = 24. 
Note that the results of Table II tell us that the es- 
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Figure 2: (Color online), (a) The entanglement entropy 
Sl(£ = L/2, 5) vs 5 for the spin-l Blume-Capel model with 
7 = 1.1 and for several values of L. (b) <1Sl(L/2, 5) /d5 vs 
5. (c) ASl(S) vs 8 for the spin-l Blume-Capel model with 
7 = 1.1. (d) Finite-size estimates of c L as function of 7 for 
three values of L. (e) Same as (b) but for 7 = 0.39. 



timated results obtained from the MVEED method are 
as good, if not better, than those coming from the CM. 
This table also includes the extrapolated results acquired 
by using the VBS method.— 

As we already mentioned, at the same time we get the 
finite-sizes estimates of 8 C L , we also obtain the finite-size 
estimates of the central charge c L . In Fig. 2(d), we depict 
these estimates for the lattice sizes L = 16, 24 and 36 for 
several values of 7(<5 C ). We clearly see, in this figure, that 
the estimates of c L tend to c = 1/2 for 7 > 0.6 and, for 
0.4 < 7 < 0.6, exhibit a quite large change with a peak in 
a value close to c L = 7/10, which is the expected value of 
the central charge for the tricritical point of the BCM. 53 

It is interesting to note that like the CM, the method 
we propose in this paper also gives us estimates of crit- 
ical quantum points, in regions where they do not exist 
(gapped regions), like the ones presented in Fig. 2(d) for 
7 < 0.4. This happens in both methods due to the fact 
that correlation length is large than the system size con- 
sidered, i. e., £ = 1/Gl > L in this region. However, 
differently from the CM the MVEED method, at least 
in this application, allows us to distinguish these fake re- 
sults. In the case of a true critical point the maximum 
of ASl, as L increases, should saturate on a finite value, 
while in the non critical case it should decrease to zero. 
In Fig. 2(e), we show ASl as a function of 6 for 7 = 0.39. 
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Table II: The finite-size estimates of 5 C L for the Blume-Capel 
model for three values of 7 obtained from the MVEED method 
and the CM. The results in parentheses are from the CM. 



For this coupling we expect a first order phase transition, 
with a small gap. It is clear that for quite small lattices 
like L = 8 and L = 12 we observe the peak of ASl 
tending towards zero. 

Our estimate for the tricritical point (7t r ,5tr) comes 
from the maximum value of ASl(j, S) (in the j-S plane). 
In Table III, we present the estimated values of the tri- 
critical couplings obtained by this procedure. We also 
show in this table the finite-size estimates of the central 
charge c L obtained with the MVEED method. Notice 
the quite good convergence towards the predicted value 
c = 7/10. One finite-size estimate of the tricritical point 
obtained from the CM is also presented. 

The evaluation of those tricritical points, by using the 
CM, is much more difficult since in this case we should 
find the crossing of two gaps in a pair of lattice sizes or 
use the crossing of the first gap in three lattice sizes. This 
means that to obtain a sequence with a given number of 
estimates we need larger lattice sizes in the CM. The 
available finite-size estimates derived from the CM (up 
to L = 8) are of the same order of precision as those ob- 
tained by the MVEED method by using the same lattice 
sizes. However, in the MVEED method we get a larger 
sequence to extrapolated, since in this case there is one 
estimate for each lattice size. For this reason, in princi- 
ple, we should obtain more precise extrapolated values. 



V. CONCLUSIONS 

We have introduced a practical and simple method, 
called maximum value of the entanglement entropy differ- 
ence (MVEED), to estimate the quantum critical points 
of quantum chains. The MVEED method explores the 



L 




fitric 


c L 


8 


0.41047 


0.91240 


0.7065 


12 


0.41405 


0.91091 


0.7025 


16 


0.41496 


0.91052 


0.7014 


24 


0.41551 


0.91028 


0.7008 


36 


0.41560 


0.91024 


0.7003 


Ref. [49J 0.41563 


0.91024 





Table III: The finite-size estimates of 7*" c and S trzc , and c L 
for the tricritical point of the Blume-Capel model obtained 
with the MVEED. The finite-size estimate of the tricritical 
couplings taken from the Ref. for sizes L = 8, 9 is also 
presented. 



distinct behaviors of the entanglement entropy at critical 
(gapless) and non critical (gapped) couplings. We made 
two applications of this method. We calculated the criti- 
cal points for the quantum Ising chain and for the quan- 
tum Hamiltonian of the Blume-Capel model. We com- 
pared, in both models, the obtained results with those 
coming from the standard crossing method (CM). Our 
results show clearly that the MVEED method gives us 
accurate estimates for small lattice sizes, similar as those 
obtained with the standard CM. However, the MVEED 
method has some advantages as compared with the CM: 

(i) Simplicity of calculation. In order to obtain a finite- 
size estimate for the critical coupling constants with the 
MVEED method we only need the eigenfunction of the 
ground state, in contrast with the CM where it is neces- 
sary to calculate the mass gaps (two different energies) 
for two distinct lattice sizes. For this reason the CM de- 
mands more computational effort, as compared with the 
MVEED method. 

(ii) Identification of the universality class. Differently 
from the CM, at the same time the critical point is ob- 
tained the central charge c, that identify the universality 
class of the critical behavior, is also calculated. In the 
usual finite-size scaling where we use the CM, the central 
charge can also be estimated by using the consequences 
of conformal invariance in the eigenspectra of the finite 
quantum chain i 54 ' 55 We need, in this case, to evaluate 
additional mass gaps, in order to estimate the sound ve- 
locity. 

(iii) Effectiveness for the localization of tricritical 
points. Due to the property (ii) the method is quite effec- 
tive to locate tricritical points since the universality class 
of critical behavior changes at the tricritical points. This 
was illustrated in the example of the BCM (section 4). In 
the usual CM the numerical effort is higher, because the 
finite-size estimates are obtained by using the two lowest 
mass gaps of two lattice sizes or the lowest mass gaps of 
three distinct lattice sizes. 

(iv) Sensibility to detect the first-order phase transi- 
tions. Gapped quantum chains with small gaps (Cl = 
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1/Gl > L), produce crossing points in the usual CM, but 
also a maximum in the proposed estimator AS^. In the 
case of the CM this can be decided only by further anal- 
ysis of the finite-size dependence of the mass gaps, and 
usually it is necessary large lattices. Surprisingly in the 
application we did for the BCM, our proposed method 
shows, already for quite small lattice sizes, a distinct be- 
havior for the gapped and non-gapped systems. Even for 
non-critical points that are close to the tricritical point 
(small mass gaps) the estimator, instead of saturating 
in a non-zero value, goes to zero, as we should expect 
in a gapped system. This is a surprise because, simi- 
larly as the CM, the proposed method is also based on 
the variations of the correlation length, which are quite 
small. However, distinct form the CM, our method shows 
a large sensitivity already in quite small lattice sizes. Fu- 
ture applications of the proposed method will confirm if 
this is a general feature of the method or an accident of 
the present application. 

It is interesting to stress that differently from several 



existing estimators, based on the entanglement proper- 
ties, the MVEED method was tested in small lattice sizes 
producing surprisingly good results. For the above rea- 
sons, we believe the use of this method may become pop- 
ular in future applications. 

In order to conclude, we mention that a generalization 
of the presented method, based on the distinct behavior 
of the entanglement entropy, to locate quantum critical 
points, can also be introduced for higher dimensions at 
least for fermionic systems with a finite Fermi surface. In 
this case, a natural estimator would be ASl/L^ 1 . 
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